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Abstract. - Nature of the magnetic hysteresis for thin films is studied by the Monte-Carlo 
simulations, ft is shown that a reconstruction of the magnetization pattern with external field 
occurs via the creation of vortex-antivortex pairs of a special kind at the boundaries of stripe 
domains, ft is demonstrated that the symmetry of order parameter is of primary importance 
for this problem, in particular, the in-plane magnetic anisotropy is necessary for the hysteresis. 



Introduction. - Traditional views assume that under the thermodynamic equilibrium 
conditions a system should either be homogeneous or consist of macroscopically large domains 
of homogeneous phases. It appears, however, that equilibrium or very long-lived metastable 
states occur frequently with a mesoscale heterogeneity (modulated phases, patterns, etc.; 
for a general review see, e.g., Ref. [1]). These inhomogeneous states may be either regular 
(e.g., stripes, stripe domains) or "chaotic". The stripes in high-temperature superconductors 
and doped Mott insulators [2,3], supramolecular self-assembly in organic chemistry [4], and 
"heterogeneous fluctuations" in metallic alloys [5] should be mentioned in this context. Stripe 
magnetic domain state of thin ferromagnetic films [6] under certain conditions turns out to 
be unstable with respect to the formation of complicated two-dimensional "chaotic" pattern 
[7-9], which provides another example of the mesoscale pattern evolution, interesting not only 
conceptually but also for the applications related to the information storage. 

It is commonly accepted now that the formation of the mesoscale heterogeneity is a result 
of frustrations in the system which can result from either geometric factors [5, 10, 11] or 
competing interactions, the long-ranged forces such as Coulomb or dipole-dipole interactions 
being of primary importance [3,7, 12-14]. General concepts of "avoided criticality" [12, 13] 
and "self-induced glassiness" [14] have been proposed to treat this situation. However, we 
still have no detailed theory to describe a formation of heterogeneous mesoscale states in 
frustrated systems. In particular, there is no satisfactory description of the experimentally 
observed [8, 9] "chaotization" of the magnetic stripe domains beyond the oversimplified Ising 
model; for the latter, several Monte-Carlo [15, 16] or phase-field [17] simulations have been 
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carried out. However, the vector character of the order parameter is of crucial importance for 
the magnetism. 

In particular, the key issues for the applications such as the pattern evolution with the 
in-plane magnetic field variation and the magnetic hysteresis cannot be formulated in the 
framework of the Ising model. These questions are investigated in the present work. 

Stripe domain formation is the result of an interplay of magnetostatic energy, energy of 
domain walls, and easy-axis magnetic anisotropy [6]. For the case of magnetic multilayers, 
the latter comes mainly from the layer interfaces. Easy magnetization axis normal to the film 
plane xy appears [9]. Further we will argue that the in-plane anisotropy is also crucial for 
the magnetic behavior of the films. To demonstrate this we consider two limiting cases, that 
are the model with the two-component magnetization vector lying in the yz-plane and the 
three-component (Heisenberg) case without any xy-anisotropy. 

Model and simulation results. - We start with the following effective Hamiltonian for the 
system under consideration: 

^ = / (i + ^ (^)' - f - h^y) + (1) 

Here m is the magnetization unit vector having either two or three components, h is an 
external field, Jx,Jy,K are exchange and anisotropy parameters, and the last term describes 
the long-range dipole-dipole interaction [7] {Q is the effective magnetic charge density and d 
is the film thickness). Note that stripes appear only if the exchange interaction is anisotropic; 
a checkerboard-like structure arises for — Jy [9]. With our choice of parameters, stripes 
are formed along y-axis, that is parallel to the external field. The domain structure arises 
for not too strong external fields since at /i — s- oo the magnetization vector is just parallel to 
the y-axis. We will consider the case of a regular stripe domain structure is stable at zero 
external field (which is observed, e.g., for the permalloy-cobalt multilayers in Ref. [9]). In 
these experiments an instability of the stripe domain structure with the external field increase 
has been observed, which can be described in terms of the appearance of disclinations at the 
domain wall lines [8,18]; further a two-dimensional "chaotic" magnetization pattern is formed. 

In our Metropolis Monte Carlo [19] simulations, a discrete-lattice (300 by 300 sites) analo- 
gous to the effective Hamiltonian iQJ is studied, in particular, the exchange terms are approx- 
imated by appropriate nearest-neighbor interactions. We present here the data for the system 
with Jx = 2^Jy = A,K — 4, = 0.1 and d ~ 10 obtained at several values of the tempera- 
ture. The parameters are chosen to produce the stripe domains with a width A of about 10 
lattice periods. In each simulation, we start with a large (saturating) negative external field 
and slowly varied it to a large positive value and then back. 

Figure n shows the magnetization curves Myih) for several values of temperature for the 
two-component (yz) magnetization. Well-pronounced hysteresis loops are clearly seen. To 
reveal the source of the hysteresis, we display snapshots of the system at several points of the 
hysteresis loop (Figure |2Jl. Inside the stripe domains, the magnetization is oriented almost 
perpendicular to the film plane, so that = 1 and niz = — 1- These areas are drawn by red 
and light-blue colors, respectively. At a domain wall, the magnetization has to pass through 
the points where niy = — 1 or my = 1. These points are shown by and dark-blue and yellow- 
green colors. Mixtures or the appropriate colors is used for intermediate orientations of the 
magnetization. 
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Fig. 1 - Magnetization curves for the system with two-component magnetization. Red, black and 
blue loops correspond to the temperature values 1//3 = 0.5, 1 and 2, respectively; other parameters 
are given in the text. Arrows indicate the critical field, h, corresponding to the nucleation of the 
vortex-antivortex pairs. The inset shows these values of h, whereas the line in the inset shows the 
estimation Ol with In = 13 



At zero field, the magnetization rotates through the negative TOy-direction for all domain 
walls. At a certain positive h, "green" parts of the domain walls with positive ruy arise. 
Further increase of the external field results in a sudden fiip of all domain walls to positive 
niy. One can see that the magnetization vector rotates over 2t: angle at the motion along a 
closed loop containing the end point of the "green" part. Thus, the magnetization flip regions 
of the domain walls are terminated by vortex-antivortex pairs. It should be stressed that here 
both the vortex and antivortex always locate at the same linear domain wall contrary to the 
Kosterlitz-Thouless situation [20] where the vortices can travel in the entire plane. 

Similar computations have been performed for the case of three-component magnetization 
(Heisenberg ferromagnet with dipole-dipole interactions). The system has a similar stripe 
domain structure. However, in the very contrast to the two-component case, no hysteresis is 
observed within the errorbar of calculation, as Figure |31 shows. 

Physical picture. - Let us first discuss the source of hysteresis for the two-component 
magnetization. Consider the energy balance of a single vortex-antivortex pair starting with 
the case h = 0. First of all, creation of the vortex itself requires a finite energy. Near the 
vortex singularity, energy is mainly contributed by the exchange (gradient) terms in Eq.(^), 
similarly to XY-model [20] . The size of this area can be estimated as a smaller value from the 
vortex-antivortex distance I and domain wall thickness a = ^ Jx/ K. A standard estimation 
with the logarithmic accuracy similar to Ref. [20] results in the expression for the energy of 
pair at /i = 0: 

E{l<a)^2TT^JJ'y\n{l/r), (2) 
E{1 > cr) = 2n^/XX,\n{a/r) 

Here the short-range cutoff of order of the lattice constant, tq, is introduced. 




Fig. 2 - Snapshots of the stripe-domain system with the two-component order parameter at several 
points of the hysteresis loop for f3 — 1. Magnetic field is h — 0, h — 0.3, and h = 0.6, from top to 
bottom shot. The inset shows the color legend for the orientation of local magnetization. 
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Fig. 3 - Magnetization curve for the case of tree-component magnetization at temperature = 1. 
All the parameters are the same as for 1//3 = 1 magnetization curve in Figure 1. 

The existence of hysteresis means that the value of E{1 > a) is large enough in comparison 
with the temperature, otherwise thermal fluctuations should destroy oriented domain walls. 
The external field makes the wall flip energetically favorable and decreases the barrier between 
the the flipped and unflipped wall. Let us estimate the height of this barrier. We will 
discuss the lower part of the magnetization curve, so that positive h tends to flip the domains. 
Consider the term —hmy as a perturbation. The zero-order approximation for m(r) can 
be constructed explicitly for I <^ a, because depolarization effects and the anisotropy term 
in 1^ can be neglected in this case. The first-order correction to E{1 <C a) appears to be 
^ (/ ''^y^^^) h-o ^ — In 2^ JxJyPh. For large I, this correction is obviously proportional to 
/, that is E{1 ^ (t) cx —Ih/a. One can see that the function E{1) at non-zero positive h has a 
maximum at ^ ^ cr or ^ ~ cr. Use the limit of small I for the logarithmic-accuracy estimation. 
We obtain that the maximum of E{1) take place at Z^^^ (x Jx/h and E{lmax) ~ T^yJ JxJy In -^rfi- 
There is a noticeable probability of the vortex-antivortex pair nucleation at the domain wall 
interval of the length L during the time T, if e-P^C-ma^) ~ Xk^ where short-time cut-off to 
is the shortest relaxation time of the system [2f]. Therefore we estimate that the vortex- 
antivortex pairs arise at 



It is worthwhile to note at this point that, strictly speaking, the applicability of the Monte 
Carlo scheme is proven only for the thermodynamic equilibrium case. However, it is clear 
that for the situation of several well-separated local minima, individual properties of those 
metastable states are also described correctly. The transition through energy barrier between 
the minima in the Metropolis Monte Carlo simulation is in fact a diffusion process with a typ- 
ical timescale to of about a single Monte Carlo sweep, therefore in our numerical calculations 
In^ « 10 15. The obtained dependence agrees well with the numerical data, as the inset 
in Fig. 1 shows. 

Possible effects of interaction between domain walls should be taken into account. To 
determine this interaction, it is sufficient to consider only the structure which is homogeneous 




(3) 
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in the y-direction. One should note first that the domain waU flip does not affect the dipole- 
dipole interaction term since this flip does not change z-component of the internal magnetic 
field. Therefore we can ignore the last term in Eq.(^ in our estimation. At zero external field, 
Eq.Q reads 

J ■ K \ 

-^V 2 ^^"^ (f>jdx — min, (4) 

where the phase (f> describing the orientation of two-component magnetization is introduced: 
m = (sin 0, cos 0). This is equivalent to the minimal action condition for a pendulum, 
if X acts as a time [22]. Qualitatively, the domain wall corresponds to a potential energy 
minimum of the pendulum (j) = 0, ±7r, whereas domains themselves are described by the 

near-separatrix motion « ±7r/2, ±37r/2, Domains with parallel and antiparallel walls 

are described by the conditions 0(0) = 0, 0(A) = and 0(0) = 0, 0(A) = tt, respectively. 
The difference between the values of the action for these two trajectories gives an effective 
interaction energy of domain walls. For the case of A/a 3> 1, the parallel wall configuration 
is energetically favorable, but the energy gain is exponentially small with a factor of e~^l'' . 
In our simulations as well as in the experiment [9] A/cr k,1 ~ 10, and this exponentially small 
term can be neglected. It would be interesting to realize experimentally a system with A « cr, 
where the interaction between domain walls would be important. 

In the case of the Heisenberg model (three-dimensional vector order parameter without xy 
anisotropy), the vortices and antivortices will not arise since the orientations of domain walls 
can be changed without a transition over an energy barrier. For example, for the domain 
profile passing from to^ = 1 to = — 1 orientation via y-direction, one can consider the 180*^ 
rotation of m in the xy-plane. The magnetization distribution m^; ~ sin 0-\/l — m^, rriy = 
cos0-\/l — 771^ interpolates smoothly between the domain walls passed through +?/ and — y 
(0 = and = tt) with the energy independent on 0, as one can see from Eq.(^. In 
agreement with these simple topological considerations, our numerical simulations do not 
show any hysteresis effects for the Heisenberg case. This demonstrates the relevance of in-plane 
magnetic anisotropy for the hysteresis in magnetic films and multilayers. This anisotropy can 
suppress a continuous rotation of the magnetization vector and therefore result in a hysteresis 
behavior. 

We investigate here a simplified model of the stripe-domain formation; in the conclusion, 
we discuss its relation with real experimental situation. Although the Hamiltonian is 
frequently used to describe the stripe domains [7] , it may by a concern that the dipole-dipole 
interactions due to the in-plane magnetization components are neglected in this approach. 
Indeed, the effective density of magnetic charges is —Q div m, and only the surface charges 
due to the discontinuity of the normal component of m at the film boundary are taken into 
account in Eq.(^. The bulk charge density is proportional to -I- which is zero inside 
the domains but not in the domain walls. Physically, the in-plane dipole moment of the domain 
wall may occur. The detailed study of these effects is out of the scope of the paper; let us 
only justify why they can be neglected in the cases under consideration. First of all, the bulk 
magnetic charges do not appear for the two-component order parameter since mx is absent and 
rUy is independent of y. For the three-component order parameter, the domain wall may have 
a dipole moment of order of Qad per unit wall length. This value should be compared with 
similar quantity of the domain, that is QKd. Therefore our consideration is valid assuming 
that cr ^ A. Qualitatively, even in this case the bulk charges introduce some anisotropy in the 
xy-plane and can therefore result in hysteresis effects for the three-component magnetization, 
but these effects should be much weaker than for the two-component case. 
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